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Abstract 

Wilson's exact renormalization group equations are derived and integrated for the rele- 
vant part of the pure Yang-Mills action. We discuss in detail how modified Slavnov-Taylor 
identities controle the breaking of BRST invariance in the presence of a finite infrared cutoff 
k through relations among different parameters in the effective action. In particular they 
imply a nonvanishing gluon mass term for nonvanishing k. The requirement of consistency 
between the renormalization group flow and the modifled Slavnov-Taylor identities allows 
to controle the self-consistency of truncations of the effective action. 
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1 Introduction 



One of the most important and most interesting problems in quantum field theory today 
is the description of non-perturbative phenomena. The main tools used so far to tackle 
these problems are Schwinger-Dyson equations and the formulation on a (finite) space- 
time lattice. The difficulties with both methods can be traced back to the fact that in 
quantum field theory effects on all length scales are involved in an essential way. 

The last years have witnessed the emergence of a new technique which may be capable 
of surmounting these difficulties. It originated much earlier in the work of Wilson et al. 

and is based on the idea that the effects of different length scales should be taken into 
account one after the other rather than all at a time. The resulting 'exact renormalization 
group equations' or 'flow equations' are formulated directly in continuous Euclidean space- 
time, thereby avoiding any lattice artifacts. Polchinski [0] was the first to present a version 
adapted to the needs of quantum field theory. His exact renormalization group equation 
determines the flow of the effective Lagrangian; by taking the Legendre transform a flow 
equation for the effective action has been obtained later 0-0. The effective action being 
the generating functional of one-particle irreducible Green functions, this formulation is 
more suited to the applications we have in mind than the older one involving the effective 
Lagrangian, which generates the truncated connected Green functions. 

Technically, the flow equation for the effective action is obtained by introducing an 
infrared momentum cutoff k which is then varied continuously. At a large scale k only 
the quantum fluctuations on scales above k are integrated out, so that the corresponding 
effective action incorporates the microscopic physics in the sense of an effective field 
theory. In the limit k —>■ the effective action becomes the usual one, generating the phys- 
ical one-particle irreducible Green functions. The flow equation now provides a differential 
equation interpolating between k = k and = 0, so given one can in principle calculate 
the physical effective action by just integrating the flow equation. Of course in practice it 
will not be possible to perform the integration exactly, since the functional F^ will have 
to be described in general by infinitely many /c-dependent parameters. What renders this 
approach useful, then, is the existence of approximations which take into account only a 
finite number of parameters and provide nevertheless a sensible description of the relevant 
physical effects. 

Historically, flow equations have been used first to simplify proofs of perturbative renor- 
malizability Later they have been employed for the description of non- 

perturbative phenomena like phase transitions, critical exponents, bound states and con- 
densates [n|]. Of course it is very desirable to apply this method to gauge theories like 
QCD. However one encounters here the obvious problem that the infrared cutoff k breaks 
the gauge or BRST invariance. Recently two different ways to cope with this problem have 
been proposed. In background fields have been introduced in order to maintain gauge 
invariance. However the resulting additional dependence of the effective action F^ on the 
background field leads to new difficulties. In [§]-[|Tni it has been suggested to use modified 
Ward or Slavnov-Taylor identities to treat the symmetry breaking induced by the infrared 
cutoff. These identities guarantee the gauge or BRST invariance of the physical effective 
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action at k = 0. We follow here the formulation presented in [|T0| 



The aim of the present paper is to demonstrate the suitability of this approach to cal- 
culations in gauge theories. For definiteness we consider SU(3)-Yang-Mills theory, but the 
concept is by no means limited to this particular case. The emphasis here is on the method, 
so we shall work with a rather simple but nevertheless non-trivial approximation for the 
fc-dependent effective action. However, the formalism is general and can be combined with 
other techniques already developed in connection with the flow equations. 

The paper is organized as follows: In section 2 we present the flow equations and mod- 
ified Slavnov-Taylor identities following but using a more concise and generalizable 



notation. In section 3 we outline how this framework can be employed to calculate the 
physical effective action, and we present the approximation we shall use in the follow- 
ing. Section 4 represents the main part of the paper and contains the results from the 
integration of the flow equations. We shall find in particular that the BRST invariance 
of the original theory provides us with the means to control the approximation we have 
introduced. Section 5 contains our conclusions and a short outlook. 



2 Flow equations and modified Slavnov-Taylor 
identities 

Throughout this paper we consider SU(3)-Yang-Mills theory in four-dimensional Eu- 
clidean space-time, although the extension to general SU(N)-Yang-Mills theories in d- 
dimensional space-time does not involve additional complications. 
We take the classical action to be 

s = I d^x { + ^d.A^^d^Al + d,e (d,d^ + gf\.Ay) 

- {d,c^ + gf\cAy) - L^\gf\cc'c^ + L^^U^A^^ (1) 

with 

= d,Al - d^Al + gf\AlAl , (2) 

where we have included the usual gauge fixing and ghost parts and coupled external sources 
to the BRST variations 

5AI = {d,c^ + gf\,Alc^)C 
dd^ = \gf\cC^cX 

5e = --d,A;C. (3) 

a 

Here is a Grassmann parameter. The invariance of S* at L = under BRST transforma- 
tions can then be expressed as 

= I d'xhA^^^ + Sc'^^ + Sc'^^] 
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(4) 

L=0 



To derive the flow equations we add an infrared momentum cutoff term 

ASk = ^A- Rk- A + c- Rk-c (5) 

to the action S. Here we have introduced a matrix notation which will be useful in the 
following. It is most easily visualized in momentum space, where the matrix product 
implies integrations over momenta as well as summations over group and Lorentz indices. 
As an example, A ■ R^ ■ A actually stands for 

Imposing energy-momentum conservation and global SU(3)-invariance we have 

RtUp,-Q) = RkAPW\27i)%P - q) 
Rt{p,-q) = Ru{pW\2t:)%p - q) . (7) 

The functions Rk^^vip) and Rkip) introduce smooth infrared cutoffs k in the gluon and 
ghost field propagators and vanish identically at = 0, 

Rk{p)^^ as . (8) 

Furthermore, 

Rk,fiu{p), Rk{p) 'rapidly' as ^ oo , (9) 

so that Rk and Rk act simultaneously as ultraviolet cutoffs in the integrals on the r.h.s. of 
(p!6|) and (pT]) below. Explicit expressions for Rk and Rk will be given later. 

The /c-dependent generating functional of connected Green functions Wk is defined by 

^w,[.j,x,x,k,l:l] ^ J T^^^^i^A, c, c)e-^^fc-^+^-^+x-+-x ^ (iq) 

where we have attached the index "reg" to the functional integration measure to indicate 
that the functional integral has been regularized in the ultraviolet. This regularization 
is necessary to give a meaning to ([T0|), and it should be independent of k and preserve 
the symmetries of S. It is not quite clear whether such a regularisation exists in a non- 
perturbative sense. In the following we shall nevertheless assume its existence in order to 
present a simple derivation of the flow equations and the modified Slavnov-Taylor identities 
(STI) and to illustrate the physical meaning of the quantities under consideration. Once 
we have obtained the flow equations and modified STI, however, we can use them without 
reference to (|1^) and thus do not rely on the existence of an appropriate regularization any 
more, as has been emphasized in In a more rigorous sense, the theory may directly 



be defined via the flow equations and modified STI 
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So let us, for the moment, use ( |T0|) in order to define the functional Wk- Since ASk 
vanishes for A; = 0, we see from ( |T0| ) that Wq is the generating functional of the physical 
connected Green functions. 

For general k we pass over to the Legendre transform with respect to the sources J, x, 

r.lA.c, -c,K,L,L] = A-J + x-c + c-x- Wt{J, x, X, K, L, L] , (11) 

where A=^, c='-^, . ^ -'-^ . (12) 

dJ dx ox 

Finally, F^ is defined by 

Ffc = ffc + ASu . (13) 

At A; = 0, f = Fq becomes the physical effective action or the generating functional of 
one-particle irreducible Green functions. 

Differentiating F^ with respect to k leads to the fiow equations. To cast them into a 
concise and generalizable form, we introduce the superfields 

$ = (c, c. A) and $ = (-c, c. A) (14) 

and (for later use) the external supersources 

Q = {L,L,K). (15) 

The fiow equations finally read 0-0 

(- 



2" 



(16) 



Here the inverse is to be taken in the matrix sense including the momenta as indices. The 
elements of the matrix of second derivatives are in our notation given by 

where i and j denote the group and (in the case of gluon fields) the Lorentz indices. TZk 
stands for the matrix of cutoff terms, 



( Rk 
^fc 
\ Rk 



written in a block matrix notation. STr denotes the supertrace, 

STr(A) = Tr(M ■ A) (19) 
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for some (super) matrix A, where 



M 




(20) 



in the above block matrix notation. Again, the trace includes summation over the mo- 
menta. 

Flow equations for the one-particle irreducible n-point functions can be derived from 
(p!6|) by taking derivatives with respect to the gluon and ghost fields as well as to the 
external sources K, L and L, setting them equal to zero after differentiation as usual. 
The complete set of flow equations is then represented by an infinite system of coupled 
differential equations for the n-point functions. It is possible to visualize them in the form 
of Feynman diagrams. The /c-derivative of a n-point function is given by a finite sum of 
one-loop graphs, where in contradistinction to the usual diagrams of perturbation theory 
the propagators and vertices are dressed. All the diagrams are infrared and ultraviolet 
finite provided the cutoff functions (^ appearing in ( pIBf ) are appropriately chosen. We 
emphasize that despite their one-loop appearance the equations are still exact. 

The introduction of the cutoff term AS^ in the action breaks the BRST invariance for 
L = explicitly. Using the invariance of the functional integration measure the symmetry 



breaking can be quantified and expressed in the form of modified STI for ||10|| : 



6Tk STk STk STk STk STi, 



SK SA 5L Sc SL 5c 



STr 



L=0 



7^i 



.(. 



6Q6<^ \6M'^ 



L=0 



(21) 

The external supersources Q have been introduced in (|15]). In ( [2T| ) the left hand side 
represents the usual variation under BRST transformations, while the right hand side 
specifies the symmetry breaking term, vanishing at A; = 0. Both sides are to be taken at 
1 = 0.0 

The modified STI ( PT| ) are similar in form to the flow equations (P^. Expanded in the 
flelds and external sources the former are also equivalent to an inflnite system of coupled 
nonlinear equations for the fc-dependent one-particle irreducible n-point functions. Again, 
the trace term can be represented as a flnite sum of one-loop diagrams with dressed vertices 
and propagators. Infrared and ultraviolet flniteness is guaranteed by appropriate choices 
for the cutoff functions. 



Another useful identity is |p|- IC 



6K 



6c 



(22) 



where d stands for differentiation with respect to the space-time variable. (ESp is obtained 
easily by performing a redeflnition of the fleld c in (p^O]). 

1. The modified STI are called "fine-tuning 



different but equivalent formulation is presented in 
conditions" there. 
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The L-dependence of f a: is trivially fixed to be 



6L 



-d-A 

a 



(23) 



The last two identities are already valid at the classical level, i.e. with replaced by 
1) then simply state the absence of quantum corrections to the classical 



S. (H) and (I 
relations when expressed in terms of 
The crucial point, as shown in UTO 



STI, as well as (^2]) and 



is that the flow equations respect the modified 
in the following sense: Once the latter are satisfied at a 
certain scale k = k, they are also satisfied for all k < k provided is obtained from 
by integration of the flow equations. In particular, at = the physical effective action 
Fq fulfills the usual STI. Hence the meaning of the modified STI for F^ is that they imply 
the BRST invariance of the physical effective action at A; = 0. 

We remark that the modified STI, as well as the flow equations, can be extended 
trivially to include matter fields. 



3 A simple approximation 

The physical quantity one is interested in is the effective action f q. Let us briefly outline 
how the flow equations together with the modified STI can be used to calculate Fq in 
principle. The general idea is to begin with a certain Fg at a scale k which has to be 
large in comparison with all physical scales of the theory, and then to integrate the flow 
equations down to /c = 0. Although the exact form of Fg is not known a priori, we 
expect from universality that the resulting Fq is widely independent of the starting point 
F^. In the present case of an asymptotically free theory the relevant parameters of F^ 
are the coefficients with non-negative mass dimension in an expansion of F^ in fields and 
momenta.0 In addition it is essential for F^ to possess the right symmetry properties. In 
our case of BRST invariance this is accomplished by imposing the modified STI. 

The main problem in practical calculations is of course the integration of the flow 
equations. Since they represent an infinite system of coupled differential equations for the 
n-point functions, an exact solution seems impossible. Instead one is forced to truncate 
the system by taking into account only a finite number of parameters representing the 
/c-dependent effective action F^. In practice one then inserts the truncated representation 
of Ffc in the r.h.s. of (0) and projects out the contributions to the flow of the parameters 
under consideration on the l.h.s. Similarly the modified STI represent an infinite system 
of coupled nonlinear equations for the n-point functions. By using the truncated form of 



Ffc in (pTf) and performing appropriate projections one arrives at a set of equations for the 
parameters representing F^ at every scale k. 

Now taking a functional F^ which fulfills the modified STI exactly and integrating the 
flow equations down to a lower scale k exactly, we would obtain some F^ which again 



^We include here the parameters which are usually termed marginal. 
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satisfies the modified STL However, due to the necessary approximations described above 
this is no longer true in practice. As a consequence we cannot expect the physical effective 
action Tq to be BRST-invariant if we just integrate the truncated flow equations. 

Let us look at this problem more closely in the special case of the gluon mass. As will 
be shown below, for a pure SU(3)-Yang-Mills theory the modified STI demand a gluon 
mass term in L^ as long as 7^ 0. Stated differently, for the physical effective action at 
A; = to be BRST-invariant, the effective action at a scale k ^ must contain a precisely 
specified gluon mass term. Due to the truncations, however, the gluon mass will deviate 
from the value prescribed by the modified STI in the course of the integration of the flow 
equations, even if one started with the correct value at a higher scale. Since the mass is 
a relevant parameter, one will thus not end up with a BRST-invariant effective action at 



as possible to exact BRST invariance within our restricted space of parameters. Since in 
an exact calculation the modified STI would be fulfilled automatically given the correct 
starting point, we can in fact turn the argument around and consider the deviations due 
to the truncations as providing an internal validity check for the approximation. 

To be more precise, we have to divide our set of parameters into 'dependent' and 
'independent' ones, where the dependent parameters are entirely fixed by the modified STI 
in terms of the independent ones. The number of independent parameters is the same as in 
perturbation theory with BRST-invariant regularization. Now we use the truncated flow 
equations to determine the flow of the independent parameters, but rely on the modified 
STI to obtain the flow of the dependent ones. On the other hand, the flow of the dependent 
parameters could also be read off directly from the flow equations without recourse to the 
modified STI. Due to the truncations, the two results obtained in this way will in general 
differ slightly. This difference can in turn be regarded as a measure of the validity of the 
approximation. 

After this general discussion let us now return to the special case of a SU(3)-Yang-Mills 
theory. The classical action including external source terms was given in (|ip. A natural 
ansatz for a truncated fc-dependent effective action would be to include all perturbatively 
relevant couplings, i.e. all couplings with non-negative mass dimension. So we choose as a 
truncation 



k = 0. 



These problems are solved if we impose the truncated form of the modified STI on L^, 
at every scale k. In this way we are 'near' to the L^ satisfying the exact modified STI 
during the whole flow and our approximation to the physical effective action Tq is as close 




(24) 



where 




(25) 
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A priori one would expect the coefficients of all the terms in (|2^) to be independent quan- 
tities. However, some of them are related, as will be explained in detail below. For similar 
reasons certain terms are omitted which would be present from dimensional reasoning 
alone. Furthermore, powers of the wave function renormalization constants and 
for the gluon and ghost fields have been extracted from the coupling constants, which will 
prove to be convenient in the course of the calculations. 

Let us first consider the coefficient of the LA-vertex, which is fixed to be the classical 
one by (^Sf ). For clarity we denote the 'classical' gauge fixing constant a by a^^^ in the 
following. It appears in the BRST transformation (Q) and through the LA-vertex enters 
in the modified STL Note that the gauge fixing constant in Tk, i.e. the coefficient of the 
longitudinal part of the gluon two-point function, will in general be fc-dependent, whereas 
Q,io) remains fixed during the flow. 

The identity (^) demands that Tk depend only on the combination dc — K, which 
implies the equality of the corresponding couplings. In particular the formally relevant 
coefficients of the ccAA- and the cccc-vertex are related to the irrelevant KcAA- and 
i^cifc-couplings respectively and are thereby also rendered irrelevant. Hence they are not 
included in our ansatz The same property of Tk prevents the appearance of a mass 

term for the ghost field.Q 

We have introduced a /c-dependent gauge fixing constant and a gluon mass m^. The 
coefficients of the AAA-, the KcA- and the Lcc-vertex are in general different and have 
been denoted by the three coupling constants g^, and gj^. We would like to remark 
that due to the relations between the different couplings implied by the modified STI the 
number of independent couplings is equal to the number of independent couplings in the 
relevant part of a BRST-invariant action. In particular we shall show in the next section 
how the modified STI fix the couplings ml, a^, g^ and gj^ in terms of g^ and the product 
Z^a^^^ for every scale k. 

As implied by (p5D, the coupling constant of the four-gluon vertex is given by (g^)"^- 
This is a further approximation which is not justified by the modified STI. In fact these 
predict a deviation of the four-gluon coupling constant from the value {g^Y- However we 
shall see in the next section that within the region of validity of our approximation the 
difference between the coupling constants g^, g^ and g^^ may be neglected for calculational 
purposes. We therefore expect the same to hold true for the deviation of the four-gluon 
coupling constant from {gkY- 

Furthermore the modified STI enforce the appearance of an additional four-gluon cou- 
pling with the group-tensorial structure 

^ab^cd + ^ac^bd + ^ad^bc ■ (26) 

The same structure is also generated by the flow equations, as it of course must be for 
consistency reasons. We neglect this coupling in our approximation because it does not 
contribute directly to the flow of the coupling constants (7^, g^ and gl^, which are our 
primary objects of interest. In fact it only contributes to the flow of the gluon mass and 

^Observe however that the cutoff function (ESh does mtroduce a mass-hke term ^ fc^ for the ghost field. 
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the four-gluon coupling itself. Concerning these additional approximations one should, 
however, always bear in mind that the comparison of the flow of the dependent parameters 
taken directly from the flow equations with the one determined via the modified STI 
provides us with a measure of the validity of the used truncation, thus indicating every 
eventual breakdown of the approximation. This concludes the discussion of the ansatz for 
Tk used in the following. 

To complete the description of the set-up used for the calculations in the next section 
we need to specify the form of the cutoff functions. These functions may, and in our case 
will, contain the fc-dependent parameters appearing in F^. We choose 




as the cutoff function for the gluons and 



Mp) = (28) 

for the ghost field. They have the demanded properties, namely to introduce infrared 
cutoffs in the respective propagators and to fulfill the conditions (||) and @, provided that 
ttfc > and limfc^o^fc ^ (typically one finds limfc^o"^! = 0, cf. section 4). 



The choice ( pTf ) may appear unnecessarily complicated at first sight. However this 



particular form is chosen to give the propagator appearing in the loop integrals generated 



by the trace terms in (|T^) and (^) a fairly simple structure. Employing our ansatz for 



we obtain for the gluon propagator 

1 f 1 _ p-ip'^+^D/k'^ ( n n \ 1 — p-(P^/"fe+'"|)/'=^ n n \ 

^ ' , U..-^ + \, ^ , Pj^Uab_ (29) 

A similar calculation for the ghost propagator yields 

1 l-e-P'/'=" 



Zk P^ 



6"' . (30) 



Due to the particular structure of the propagators we are able to perform all loop integrals 
appearing in our approximation analytically. However, elegance often comes at a price. In 
our case the drawback is that the cutoff terms involve the fc-dependent parameters Z^, 
Z^, m| and a^. Hence on the r.h.s. of the flow equations the ^-derivatives of all these 
parameters appear on evaluating dklZk, as will be seen explicitly in the next section. 

Note that the functions appearing in (^) and ( |30D as well as the coefficients of the 
transverse and longitudinal momentum projectors in (|27|) and (p9D are analytical functions 
of p^, even if m| < 0. In particular they do not contain poles. The limits — > and 
ml ^ in (pTl) are, despite appearance, well-defined and interchangeable. This guarantees 
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that the functional = f ^ + ASk is local for all values of m^, i.e. the n-point functions 
can be expanded in powers of the momenta, which should be the case for a sensible infrared 
cutoff. 

One might be tempted to choose different mass parameters in the transverse and the 
longitudinal parts in ( p4D and ( ]27|) and hence also in (pQ]), since it is only the longitudinal 
part which will be fixed by the modified STL However, in this case the coefficient of the 
non-local term p^Pu/p^ in O) and (pT]) would not vanish in the limit p^ — > 0, and 
would not be local. 



4 Solving the truncated flow equations 

We shall now turn to the integration of the flow equations using the concepts and the 
approximations explained in detail in the previous section. Employing the ansatz ( p4D the 
flow equations and the modified STI take on a characteristic form which we shall sketch 
in the following. 

We expand the modified STI in the fields and external sources to obtain equations for 
the one-particle irreducible n-point functions. Expanding these in turn in powers of mo- 
menta and performing appropriate projections we arrive at the following set of equations: 

/f (/ifc, Oik) + ^^^(0^ /2' (/^fc, afc)^ 

/l^(/ifc, ttfc) + ^r^(o) /2^(/^fc, j 

fl{^J'k, ttfc) + 2A^(0) f2if^k, ttfc)^ • (31) 

Here we have introduced the dimensionless mass parameter fik = fnl/^"^- The /j are 
complicated functions of fik and and do not depend explicitly on k. Hence k enters in 
the system explicitly only through the relation between ml and /z^. 

Furthermore we have made still another approximation by equating the three coupling 
constants g^, and g^ on the r.h.s. and denoting them collectively by gk- Without this 
approximation we would have to replace the r.h.s. of (^) by a sum of similar expressions 
with gl replaced in turn by the products {gkY, gk9k ^ ... of the three coupling constants. 
The factor of three in front of gl reflects the number of colors, while l/(47r)^ originates 
from the integration measure of the loop integrals. 

We see that in the perturbative limit of vanishing coupling constants gk ^ ^ the 
modified STI lead to the usual relations m| = 0, Z^joik = and gk = gk = 9k- The 

last two equations in (pi]) with the three coupling constants set equal to gk on the r.h.s. 
can thus be considered as the first step in an iterational procedure to determine the values 
of g-^, g^ and g^. What is more, these two equations make it easy to check whether the 







^' = ^ = 


{Any 








(47r)2 


9k -9k 


^9l 


9^ 


(47r)2 


9k -9k 


^9l 


9^ 


(47r)2 
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approximation of equating the coupling constants is admissible. In fact we shall find this 
approximation to be valid in the region of parameter space we are interested in, which in 
turn is limited by the validity of the general approximation (|2^) . Since we consider an 
asymptotically free theory and the starting scale k is supposed to be large as compared 
with the physical scales of the theory, e.g. the confinement scale, we have (yf^ <C 1 and the 
approximation of equal coupling constants is certainly acceptable for the starting point T^. 

In the same way we could have treated the four-gluon coupling constant as different 
from (g^)"^ in the ansatz and then relate it to {g^y through an equation similar to 



(pip. However, we expect the difference between these quantities to be negligible within 
the region of validity of our approximation. 

We observe that enters in (|31|) only through the combination Z^a^^\ and does 
not appear at all. This simplification was in fact the reason for extracting powers of 
the wave function renormalization constants from the couplings in p^). By counting the 



parameters and the equations — regarding g^, g^ and gl^ as different for the moment 
— we find that the modified STI in the form (^Tj) fix /i^, a^, g^ and g^ in terms of g^ 
and Z^a^^\ Thus apart from the wave function renormalization constants we end up 
with two independent parameters, corresponding to g and a in perturbation theory with 
BRST-invariant regularization. 

From now on we shall always use the approximation of equal coupling constants, i.e. 
9k = gt — 9k — 9k unless otherwise stated, so we only need to consider the first two 
equations in (0). Given and we obtain a system of quadratic equations for g1 
and 1/Z^a^^\ which can be solved easily. On the other hand we can calculate Hk and 
numerically as functions of gt and Z^a^'^\ 

Before we come to the discussion of the solutions of ( PT| ) let us consider the limit 
(^(0) 0, i.e. the Landau gauge. One can show that in this limit the second equation in 
(0) implies ak/Z^a^'^'^ = 1 irrespective of the values of gk and /i^. Hence = and the 
first equation in (|3lD becomes 

^^k = -^,{fi'{f^k,0) + mf^,M , (32) 

which is much easier to work with. Also the functions /j take on simpler, though still not 
enlightening forms in this limit. 

We shall now turn to the discussion of the solutions of (^1]). It is easy to obtain solutions 
for small gk iteratively. Neglecting terms of the order of g^, the first two equations in (^) 
become 



(47r)2 8(a, + l) 
1 1 _ 3gl 7- lOak - hal 

'^k~Z^ ~ (47r)2 48(afc + l)2 



(33) 



For brevity we have used the quantity on the r.h.s. However it can be replaced by 
in the given order in g^ as follows from the second equation in (|33D . 
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It appears problematic at first sight that the mass parameter /i^ becomes negative for 
ctfc < 1. However in the hmit of vanishing momentum the gluonic two-point function 
obtained from becomes 

which is well-defined and non-negative for all values of fik- Since the propagator ( ^9]) 



is given by the inverse of this two-point function, the effective gluon mass is seen to be 
non-negative. 

We shall refer to ( ^3l) as the one-loop STI in the following since they are valid up to 
second order in gt- We remark that our expression for fik differs from the one-loop result 



obtained in [|10[, where a different cutoff function has been used. Hence the one-loop STI 
are non-universal in the sense that they depend on the choice of the cutoff terms. 

The solutions of (|31|) for general values of Qk are shown in figure 1, where we have 
plotted fik as a function of the strong fine-structure constant as = fi'fc/47r for fixed values 
of ak- The reason for using ak as a parameter instead of Z^a^^^ is just technical, but 
since both quantities are /c-dependent, the use of Z^a^^^ would not be more advantageous 
anyway. The values of ak depicted in figure 1 are, from bottom to top, ak = 0, 1 and 3. 

We shall see below that our approximation breaks down at as ~ 1.4, so effects showing 
up at much greater values of as are considered as potential artifacts of the truncation. We 
finally remark that /ifc does not vanish for gk > in the 'Feynman gauge' = 1, although 
it does so in the one-loop approximation (|33D. 

After the discussion of the modified STI let us now turn to the fiow equations. Expand- 
ing (|TB|) in the fields and external sources we obtain flow equations for the one-particle 



irreducible n-point functions. Some of these are represented diagrammatically in figure 
2, where our ansatz for Tk has already been taken into account in that on the r.h.s. all 
vertices not appearing in (p^) are omitted. Since the cutoff functions depend on k not only 
explicitly, but also implicitly through the parameters Z^, Z^, ml = fi^k'^ and ak, we have 

O T-> dRk uu , nu o 1 ryA , ^ Rk nu o , dRk uu 

dkRk,t,u = r,j + ^, lyA Ok In Z^ + ^ dkfik + Okak 
ok alnzf oiik oak 

on the r.h.s. of (0). Now we employ our ansatz for tk, expand the flow equations for the 
n-point functions in powers of external momenta and perform appropriate projections to 
arrive at equations for the parameters of F^. Using the approximation of equal coupling 
constants Qk = Qk = 9k = 9k discussed above, the flow equations for the 'independent' 
parameters Z^, Z^ and gk flnally take on the following form: 

dlnXk 39k fr\f ^ i\f .d\nZ^ .dlnZu 

ah^ = (i^(^o(/^-«^) + ^>-"^)^ + ^^(^-«^)^ 

+ h^^{fik,ak)-^^ + h^{fik,ak)-^^] , (36) 
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where represents the independent parameter. For definiteness we choose to take the 
equation for din g^/ din from the flow equation for the ccA- or equivalently the KcA- 
vertex. The functions /i* appearing above depend only on the parameters /i^ and and 
not explicitly on k. The terms involving din Z^/ din k^ etc. originate from the expressions 
for dkRk,^iv and dkRk in (|35D. Observe that the classical gauge fixing constant does 
not enter in fR3). 



On the r.h.s. of (|36D the fc-derivatives of /i^ and ak appear. Since these parameters are 
given in terms of Qk and Z^a*^''-* through the modified STI, one can obtain expressions for 
d fik /d In k"^ and dak/ din k^ by taking the fc-derivative of the first two equations in ([31|). 
We write them in the form 



din k"^ 



dlngk 3gl ( ^ r^-. dlnZ^ 

+ h>;,{f,k, ak, + h^^if^k, ak, Zta^'^) 



91nP ^ '51nfc2 



dak „ <91nZ^ , ^/ (91n^fc , 3^^ ^o^.SlnZ^ 



ainfc2 "91nP V ^ /^InP (47r)2 \ ^ ' dlnk'^ 

+ /.^(/..,a.,Z,V°))^ + /.:(^„a„Z,^a(o))^) . (37) 

The functions /i*, now depending additionally on Z^a^^\ are given by the functions /j 
appearing in (^), their derivatives with respect to yU^ and a^, and some factors of /Xfc, «fc 
and Z^a(°). 



and ( P7| ) together represent a system of five linear equations for the fc-derivatives 
of the five parameters Z^, Z^, gk, iik and a^, where the coefficients are functions of gk, 
Hk, ak and Z^a^^\ Since only two of the latter parameters are independent by virtue of 
the modified STI, we conclude that the fc-derivatives of all parameters are fixed by two 
independent quantities, say gk and Z^a^^\ 

Choosing initial values for gi, a^^"^ (and Z^ = 1, Z| = 1), and determining the 'de- 
pendent' parameters /i^ and aj^ with the help of the modified STI, we can now integrate 
the fiow equations using ( P^D and (^7\). We remark that the division into dependent and 
independent parameters is not unique and that we could as well regard fik and ak as 
independent parameters determining the values of gk and Z^a^^\ 

Before we present the results of the integration of the flow equations let us, however, 
have a look at the way we can control the consistency of our approximation. In (|37D we 
have taken d^k/ dlnk"^ and dak /d In k"^ from the /c-derivative of the modified STI, but 
alternatively we could use fiow equations derived from (0) for dfik/dlnk'^ and dak /d In k'^ 
(as in the case of the other parameters). Due to the used truncations the corresponding 
results will generally be different, and this difference can be taken as a measure of the 
validity of the approximation. 

Concerning the additional approximation of equal coupling constants, we could have 
taken the equation for dlngk/dlnk"^ in (^) just as well from the flow equation for the 
three-gluon or the Lcc-vertex. If the approximation g^ = 9k — 9k hold, we should 
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obtain at least approximately the same result for din gk/ din k"^ irrespective of the equation 
we use. This represents an additional condition to be fulfilled within the region of validity 
of our approximation. 

Let us briefly discuss the Landau gauge a^^^ = 0. In this case is fixed to be zero 
for all values of k by the modified STI, as was remarked above. Hence the modified STI 
imply that = is a fixed point in the Landau gauge. This is consistent with the flow 
equations, which demand that the /c-derivative of ak vanish at ak = 0. 

We now turn to the solutions of the flow equations and start with the discussion of the 
one-loop approximation, i.e. we neglect terms of the order of gf.. Again, it is easy to solve 



and (|37D iteratively in the limit of small gk, and we find 

dlnZ^ _ 3gl 13 - Sa^ 

dlnk^ ~ (47r)2 6 

glnZg ^ 3gl 3-ak 

dlnk^ (47r)2 4 

d In gk ^ 3gl 11 
ainP (4^)2 g 

dfik 



dink 



2 







dak dlnZ^ 



"^ITTTTT , (38) 



dlnk"^ dlnk^ 

where we have made use of the one-loop STI (^) in order to replace fik and Z^a^^^ on 
the r.h.s. The first three equations represent the results familiar from perturbation theory, 
while the last two state that the modification of the STI does not affect the /3-functions 
at the one-loop level (remember that ak = Z^a^^^ in perturbation theory with BRST- 
invariant regularization) . 

Concerning the consistency of our approximation, we find that the flow equations for 
dfik/d In k'^ and dak/ din k"^ derived from (|T^) reproduce, to second order in gk, exactly the 
results in (0). Here we have used (^) in order to rewrite the r.h. sides. Hence the flow 
equations are entirely consistent with the modified STI at the one-loop level. Likewise 
the expressions for din gk/d In k^ taken from the flow equations for the three-gluon and 
the Lcc- vertex are identical with the one in ( ^Sf ) in second order in gk- This is of course 
necessary for our approximation to be valid within the range of parameters where one 
expects the one-loop approximation to be reliable. 

Let us now go beyond the one-loop level. Before we present the full solution, however, 
we would like to discuss briefly a different approximation that makes an analytical solution 
of the flow equations possible. It consists in arbitrarily setting fik equal to zero for all 
k, thereby ignoring the modified STI. Comparing with the full solution below where the 
modified STI are taken into account, we can thus estimate the effect of neglecting the gluon 
mass parameter. Furthermore we work in the Landau gauge. 
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Setting fik and ak to zero, ( ^6|) becomes 

d\nZ^ _ 3gl /13 l + 1501n(4/3) d\nZ^ l-361n(4/3) d\nZ^\ 

din P ~ (47r)2 \ 6 36 din F 216 Slnfc^ J 

glnZg _ 3gl /3 3 ln(4/3) (91nZ^ 31n(4/3) dlnZ^ \ 

91nP ~ (47r)2 (^4 ^ 4 SlnP ^ 4 Sin P j 

ainc/fc _ 3gl ( 11 l + 2041n(4/3) dlnZ^ 1 - 3601n(4/3) dlnZ^ 

dlnk^ ~ (47r)2 [~~Q 72 SlnP ^ 432 din k\ 



(39) 



In second order in we recover the one-loop results ( |3q) with = 0. 
From the solution of the system of linear equations ( pOf ) we obtain 



dlngk 2.6 •3(7i/(47r)2. (2.7 -3^7^/(4 



IT 



dlnk^ (4.8 - 3^7,7 (47r)2) (0.81 - 3gl/{4ny) 



(40) 



where we have given approximate numerical values for the appearing numbers. 

Since SU(3)-Yang-Mills theory is asymptotically free, we start with a small coupling 
constant gj;. at some large scale k. Lowering the scale the coupling constant increases 
according to the solution of (^0]) until it reaches the point Sgl/^iny = 0.81, where the 
/3-function given by diverges. At this point fc^ as a function of gl reaches a maximum, 
which means that we cannot integrate the equation beyond this point. In contrast to 
perturbation theory, where we run into a Landau pole, the solution here just ends at a 
finite value of gf. 

In figure 3 we have plotted as = gl/^'^ as a function of In(fc2/A;2) in this approximation 
(dot-dashed curve), together with the three-loop perturbative result (dotted curve) and 
the full solution to be discussed in the following. 

Taking the modified STI into account, we determine the flow in the two-parameter 
space {{ak — 1)/(«a: + l))5'fc/4vr) (which maps values of ak between zero and infinity into 
a finite interval). As discussed above, the fc-derivatives of all parameters are already fixed 
by these two independent ones. One could have chosen another set as well, e.g. Z^a^^^ 
instead of the above function of ak, but the present choice turned out to be particularly 
convenient. The flow towards lower values of k in this parameter space induced by the 
flow equations is indicated in figure 4 by a vector field, where each vector represents the 
negative of the derivative with respect to In k"^. To make the directions and the magnitudes 
of all the vectors in the diagram visible, we had to rescale the vectors in a slightly unusual 
way, where the actual lengths L are replaced by 0.06 + 0.015 • L. 

Starting with small values of gk because of asymptotic freedom, we see from figure 4 
that the flow is directed towards larger values of gk, as expected from perturbation theory. 
The flow in the horizontal direction depends on the initial value of ak- For ak < 13/3 
the system is driven towards smaller values of ak, for ak > 13/3 towards larger ones. The 
specific value 13/3 can be read off from the one-loop formula (pSf). As shown in figure 4, 
this behaviour remains unchanged for larger values of gk- 

The curves appearing in the figure represent the boundary of the region where our 
approximation is valid. The different curves are obtained by imposing different constraints 
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on the deviations resulting from our approximation as discussed above. The approximation 
can be considered vahd below both curves and unreliable above. More details on the 
constraints determining the curves will be given below. 

We find that for initial values > 13/3 the system fiows out of the region of validity of 
the approximation quickly, the values of Qk still being small. Therefore we concentrate in 
the following on initial values ak < 13/3. Here the fiow converges asymptotically towards 
ak = 0, the Landau gauge case. So = is not only a fixed point of the fiow, but even an 
attractive one. This is a very important feature because it shows how independence of the 
gauge parameter emerges in our formalism: Starting with an arbitrary value < 13/3, 
corresponding to a certain value of Z^a^'^\ the system fiows towards ak = 0, and hence at 
a lower scale k the functions /j and /i* appearing on the r.h.s. of the modified STI and the 
fiow equations approach the form they have in the limit — 0. Thus in the limit A; — i- 
the parameters in the effective action become independent of a^^\ provided this feature 
persists beyond our reliability bounds. 

Let us now discuss the curves in figure 4 representing the boundary of the region of 
validity of the approximation. From the above discussion it is clear that this region is 
determined by several conditions. First of all the deviations of the coupling constants 
and gj^ from g^ must be small. We can estimate the relative deviations from the last 
two equations in (0). Furthermore the fc-derivatives of the parameters must not change 
drastically when we calculate them in one of the alternative ways described above. We 
concentrate here on the /9-function of the coupling constant, i.e. on din g^/ din k"^. We 
require its relative deviation to be small when we use the fiow equation for the three-gluon 
or the Lcc-vertex instead of the one for the ccA-vertex in (|36D , and likewise when we take 
dfik/dlnk"^ or dok/dlnk^ from (|TB|) instead of using the /c-derivative of the modified STI 
in (^TD . Now the solid curve in figure 4 represents the boundary of the region where all the 
relative deviations mentioned above are smaller than 3%, while for the region bounded by 
the dot-dashed curve they are smaller than 10%. Both curves consist of several smooth 
parts, each of which represents one specific constraint. 

The region bounded by the dot-dashed curve extends to gl/4:TT ^ 2.8 at (a^ — l)/{ak + 
1) ^ 0.2. However to access this region starting with a small value g^ one would have to 
fine-tune the initial value a^. The most interesting region due to the attractivity of the 
fixed point = is the one with small values of ak- Here the strongest constraints are 
obtained from replacing the equation for din gk/ din k^ in ( PB]) by the one taken from the 
fiow equation for the three-gluon vertex for the solid curve and from replacing the equation 
for d fik /d In k"^ in ( P7D by the one taken from (|TBp for the dot-dashed curve. Hence the 
point where our approximation is expected to become inappropriate is determined by the 
value of gk where the fiow equations and the /c-derivative of the modified STI cease to be 
consistent. Figure 5 shows the /^-function of gk at = as calculated from (|36D and ( pTj) 
(dot-dashed curve) and by taking dfik/dlnk"^ from ([T6|) instead of using the ^-derivative 
of the modified STI (solid curve). More precisely, dlnas/dlnk^ is plotted as a function 
of as, where as = gl/^'^- 

The result of the integrated fiow equations is shown in figure 3 in the case of the Landau 
gauge ak = 0. Since ak = is an attractive fixed point, the curves look the same for all 
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initial values a^, provided the starting scale k is large enough. In figure 3 we have plotted 
gl/4:7r as a function of ln(A;^/fc^), so the fiow is directed towards the right. The initial value 
for the coupling constant is gj^ = 1. The figure shows, from left to right, the perturbative 
three-loop result |T^, the full solution from the fiow equations, and the analytical solution 
in the approximation /i^ = 0. We see that the three curves are nearly identical as long 
as as = 5'fe/47r is smaller than 0.2 approximately, but already for 0:5 ~ 0.5, where our 
approximation is still expected to be reliable, deviations become visible. We conclude that 
the effect of the non-vanishing gluon mass parameter becomes important here. Observe 
that the value of k"^ where our approximation ceases to be reliable, i.e. where as ~ 1.4, 
lies beyond the Landau pole of three-loop perturbation theory. 

Finally, let us have a look at the behaviour of the gluon mass parameter, again for 
ak = 0. Figure 6 shows the normalized mass parameter m\/k'^ as a function of ln(A;^/A;^). 
It is driven to zero during the fiow, and we emphasize again that this is nontrivial, since a 
slightly different starting point mi is not expected to lead to a similar behaviour. 



5 Conclusions and Outlook 

To summarize, we have shown that the method proposed in [1^ to use flow equations for 



the description of gauge theories is well suited for practical calculations. This has been 
demonstrated here for the case of a SU(3)-Yang-Mills theory in a simple but non-trivial 
approximation, where the fc-dependent effective action is approximated by its perturba- 
tively relevant part. 

We have emphasized the crucial role played by the modifled STI in this approach both 
for the general concept and for practical calculations, where approximations are inevitable. 
These identities allow to control the used approximation by providing internal consistency 
checks. We believe that this is a new and very useful property, which of course rests upon 
the local gauge symmetry of the original theory. 

We have reproduced the usual perturbative one-loop results in the limit of small cou- 
pling constants, and we have studied in detail the behaviour and the effects of the gluon 
mass term (in particular on the higher order results), which is required by the modifled 
STI. Finally we have argued how independence of the gauge flxing constant arises in our 
formalism. 

The fact that the internal consistency checks did not allow us to go much beyond 
standard perturbation theory is an obvious concequence of our truncation of the effective 
action, not of the method per se: The discrepancies between flow equations and modifled 
STI for small k tell us that higher dimensional operators are required in order to obtain 
a consistent description in this regime; this is not astonishing in view of the expected 
dynamical effects such as glueballs and condensates. Higher dimensional operators should 
allow, e.g., to describe an arbitrary momentum dependence of propagators and vertices. In 
the case of non-gauge theories, such more general parametrizations of effective actions have 
already been used for studies of phenomena like the formation of bound states, condensates, 
and other non-perturbative effects [Tl|. The obvious next steps are now to investigate 
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these phenomena within non-abehan gauge theories, using generahsations of the method 
developed in the present paper. 
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Figure captions 

Figure 1: The dimensionless mass parameter fik = ^\/k'^ as a function of the strong 
fine-structure constant as = g'il'^T^ foi' different values of the fc-dependent gauge 
fixing constant a^- The values depicted in the figure are = (solid curve), 1 
(dot-dashed curve) and 3 (dotted curve). 

Figure 2: Flow equations for the one-particle irreducible n-point functions. The vertex 
functions are to be taken from F^, while the propagators are given in (pQ]) and ([30|). 
The figure shows the flow equations for the gluon two-point function, the ghost 
two-point function and the ccA-vertex. 

Figure 3: The strong fine-structure constant as = fl'fc/47r as a function of the scale k 
in the Landau gauge = 0. The initial value sX k = k is taken io he gj, = 1. 
The figure depicts the perturbative three-loop result (dotted curve), the full solution 
from the flow equations (solid curve) and the solution in the approximation /i^ = 
(dot-dashed curve). 

Figure 4: Representation of the flow in the two-parameter space ((a^ — 1)/ (a/c + l),as = 
gl/Air). The arrows are given by the negative of the derivatives with respect to In /c^, 
their lengths L rescaled to 0.06 + 0.015 • L in the figure. The curves represent the 
boundary of the region of validity of the approximation. The solid curve results from 
the 3%-constraints and the dot-dashed curve from the 10%-constraints (see section 
4 for details). 

Figure 5: The /3-function (3 [as) = d\nas/d\nk^ as a function of the strong fine-struc- 
ture constant 0:5 = g'i/^'n in the Landau gauge ak = 0. The dot-dashed curve 
results from (^) and (pTf), while for the solid curve d^k/dlnk^ is taken from ([T6|). 
The deviation can be used to determine the region of validity of the approximation. 

Figure 6: The normalized mass parameter ml/k^ as a function of the scale k in the 
Landau gauge a^ = 0. As in figure 3, the initial value ai k = k is determined by 
9-k = 1- 
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